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.•BSTRACT 

licistiiie  navigation  techniques  usins  satellite  do.;ijler  shift 
meaoureinents  are  sensitive  to  observer  velocity  errors..  model  is 
developed  for  determining  the  position  and  velocity  of  the  observer  with 
the  result  that  the  na.ig  n  process  will  be  independent  of  external 
velocity  source  error,  liie  position  and  velocity  determination  will  be 
feasible  provided  that  precision  reference  oscillators  with  instabilities 
of  about  2  p;irts  in  10^^  be  used  in  the  satellite  and  the  navigator's 
equipment,  Relativi;-'tic  and  tropospheric  effects  must  be  accounted  for 
in  the  model.  The  feasibility  also  depends  on  the  satellite  velocity 
error  being  less  than  0.05  meter  per  second. 
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I.  INTHODUCTION 

Existing  navigation  techniques  using  the  measiu:ed  doppler  shift 
of  a  satellite  signal  with  respect  to  an  observer  moving  on  the  surface  of 
the  earth  or  at  a  measured  altitude  require  initial  estimates  of  observer 
position  and  velocity.  In  these  techniques  the  position  of  the  navigator 
is  computed  by  least  squares  estimation  at  some  epoch  and  requires  the  use 
of  velocity  estimates  from  other  sources  which  are  assumed  to  be  correct. 

This  paper  demonstrates  the  feasibility  of  determining  both  position  and 
velocity  from  the  satellite  doppler  shift  data. 

The  doppler  shift  is  measured  simultaneously  on  two  coherent, 
harmonically  related  frequencies.  A  continuous  real  time  correction  fcjr 
ionospheric  refraction  is  made  to  first  order  [!]>  [2],  [3]  and  [1^].  The 
satellite  borne  oscillator  also  serves  as  the  master  clock  for  the  navigation 
system.  A  precise  navigation  program  requires  that  the  navigation  system 
clock  be  synchronized  with  the  satellite  clock.  System  concepts  of  satellite 
doppler  navigation  are  given  in  considerable  detail  by  Guier  and  Weiffenbach 
[1],  [2],  [3]  and  [It-]. 


Navigation  programs  in  current  use  determine  only  two  state 
variahlfiB,  latitude  and  longitude,  and  an  additional  parameter,  frequency 
bias.  The  frequency  bias,  f^,  is  the  difference  in  the  frequencies  of  the 
satellite  oscillator  and  the  receiver  reference  oscillator.  Typical  f 

g 

detennlnations  in  a  navigation  fix  range  from  one  part,  in  10  to  one  part 
in  10^.  The  f^  terms  must  be  removed,  or  else  be  precisely  known,  during 
the  navigation  process  or  an  error  will  be  induced  in  the  latitude.  The 
stability  pf  existing  oscillators  is  such  that  the  f^  term  can  be  represented 
as  a  constant  during  the  time  of  the  satellite  peiss.  That  f^  may  be 


CONFIDENTIAL 

doounrat  cmtaim  iMoriMtion  «lhclii^  the  nelionel  defenw  the  United  Stetee  within  the  nweninf  ol  the  Eipioneie  Lewe.  Title  18 
UeS.C.,  Section*  T93  m.  Tkc  tranemiieion  or  the  revdetion  ol  Hi  conltnli  in  any  maiiner  to  an  unaalhoriaid  penon  it  prohtbilod  by  law. 


-2- 


^  TMC  JOHNS  HOMINS  UNIVtSSITY 

APPLIED  PHYSICS  LABORATORY 

SH.Vt«  SfWMa.  HAMVLAtm 


CONFIDENTIAL 


considered  constant  implies  that  the  doppler  residual  due  to  this  term  is 

an  even  function  with  respect  to  the  time  of  satellite  closest  approach, 

It  is  shown  on  Figure  1  that  the  doppler  residual  due  to  a  latitude  error 

is  also  an  even  function  with  respect  to  t  .  The  f  terms  must  be  removed 

c  7 

or  the  fitting  process  will  misconstrue  the  frequency  bias  as  a  latitude 
error. 

The  major  \mcorrected  error  source,  when  it  occurs,  in  non¬ 
stationary  navigation  is  due  to  the  use  of  inaccurate  velocity  data  in  the 
computing  procedure.  The  effect  of  input  velocity  errors  in  current 
computational  methods  is  clearly  shown  on  Figures  2,  3,  and  5.  The 
magnitude  of  the  position  error  is  dependent  on  the  magnitude  of  the  velocity 
error  and  on  the  relative  geometry  of  the  navigator  observing  the  satellite. 
The  sign  of  the  position  error  is  dependent  on  observer-satellite  geometry 
and  the  direction  of  satellite  motion  relative  to  the  observer. 

Figures  2  and  3  show,  res:^ctively,  the  magnitude  of  longitude 
error,  b\  in  nautical  miles,  due  to  an  input  velocity  north  error,  in 
knots,  and  the  second  order  latitude  error,  &<?  in  nautical  miles,  due  to 
&Vjj.  Tlie  abscissa  of  Figures  2,  3»  ^  and  5  is  given  in  terms  of  the 
elevation  angle  at  time  of  satellite  closest  approach,  measured  with  respect 
to  the  topocentric  horizon.  Figures  4  and  5  show  the  magnitude  of  latitude 
error,  5cp,  due  to  an  input  velocity  east  error,  -bV,,  in  knots,,  and  the 
magnitude  of  the  second  order  longitude  error,  6X,  due  to  &V_.  The  position 

£1 

error  curves  due  to  velocity  error  represented  by  Figures  2,  3,  4  and  5  were 
obtained  in  a  computer  simulation  and  have  been  -verified  by  analytic  methods. 
The  position  errors  due  to  velocity  errors  are  linear  over  a  large  range, 
about  20  knots,  of  velocity  errors.  The  modeling  necessary  to, determine 
position  and  velocity  from  the  satellite  doppler  data  will  now  be  outlined. 
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Fig.  2  Longitude  Error  per  Velocity  North  Error 
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II.  DISCUSSION 

The  expression  for  the  doppler  shift  of  the  satellite  signal  as  seen 

by  the  observer,  including  general  relativistic  effects  [5],  [6],  [7],  thru 
2 

order  c  is  given  by  (l) 

(1)  ^(t,T)  = |p(t,T)  •  p(t^T)  +~j^(p(t,T)-  Rg(t))^+{p(t,T)-R^(t))(p(t,T)*rjj{T)) 

where  t  =  satellite  clock  time 

T  s  observer's  clock  time 
=  satellite  position 
=  observer's  position 
pCt,-)  =  Rg(t)  - 
p(t,T)  s  Ht,T)/lp(t,T)l 

c  =  velocity  of  light 
G  =  Newtonian  gravitation  constant 
Mg  =  earth  mass 
f^  ■  source  frequency 
Af  =  doppler  frequency 
T  =  t  +  |p(t,T)l/c 

and  the  dot  (')  indicates  time  derivative. 

Equation  (l)  is  the  modeled,  or  theoretical,  doppler  frequency 
to  be  computed  in  the  navigation  process.  Relativistic  corrections  are 
not  necessary  for  latitude  and  longitude  position  determinations.  The 
relativistic  terms  must  be  included  for  position  and  velocity  fixes,  since 
velocity  determination  from  doppler  data  requires  that  the  frequency  be 
relatively  bias  free,  as  will  be  shown  later.  Examination  of  the  terras 
of  order  c  in  (l)  reveals  that  these  terras  contribute  a  frequency  bias 
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9 

term  of  the  order  of  one  part  in  10^  for  an  orbit  with  a  semi -major  axis  of 
about  l.l6  earth  radii.  The  maximum  doppler  shift  relative  to  the  source 
frequency  is  about  2  parts  in  1(P. 

The  first  order  approximation  of  (l)  is 

f 

(2)  Af(T)  =  -  P  *  (Rg  “ 

It  is  clear  from  (2)  that  errors  in  the  observer's  velocity  and 

position  contribute  to  errors  in  the  doppler  frequency.  It  is  also  evident 

from  (2)  that  the  satellite  ephemeris,  R  and  R  is  used  as  the  reference 

s  s 

in  the  determination  of  the  doppler  frequency.  Any  errors  in  the  ephemeris 
will  contribute  errors  to  the  observer's  position  and  velocity  determination. 

The  least  squares  adjustment  of  the  observer's  position  and  velocity 
to  measured  doppler  points  n  and  other  apriori  data  is  formulated  by  (3). 


(3) 


N 


F 


-if* 

~  U 


(1=1 


t£{t  ,  T  lx)  -  Af  (t  ) 
(i^  (i'  V  |a 


where  Af(t  ,  t  |x)  represents  the  theoretical  doppler  data,  Af  (t  )  the 

(^  (1  V  [J« 

measured  doppler  data  and  x  the  state  vectors  to  be  minimized.  The 
components  of  the  state  vector  eire 


x^  =  cp,  observer's  latitude 

x,^  =  observer's  longitude 

d 

x^  =  9,  observer's  angular  velocity  along  a  meridian 

« 

X|^  =  X,  observer's  angular  velocity  along  a  circle  of  latitude. 

The  remaining  components  of  the  state  vector  r^j  and  f'jjj  will  be  constraints 

on  the  above  formulation;  details  for  computation  will  be  given  later. 
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Expressing  the  theoretical  doppler  shift  by  a  Taylor  series  to 
linear  terms  we  have 


4  5Af(t  ,  T  lx  ) 

(I.)  tJx)  =  ai. 


where  the  subscript  e  represents  the  initial  estimates.  Define,  for 
convenience 


SAf(t  ,  T  X  ) 
,  _  '  u  e'' 

b .  = - ‘-T — - 

1  ox. 


and 


R  s  Af(t^,  doppler  residuals. 

The  least  squares  adjustment  may  now  be  written  as 


The  usual  minimization  of  F  by  =0,  i  =  1,  2,  3>  ^  results  in  the 

Ji 

normal  equations  (6). 


(6) 


N 


Rb, 


N  4 

Z  E 

p=l  i=l 


b. 

1 


Ax. 

1 


=  0, 


^  =  1,  2,  3,  4 


Equation  (6)  will  be  solved  for  Ax^.  An  improved  estimate  of  x^ 
thereby  an  improved  estimate  of  the  theoretical  doppler  value  can  be 
computed.  The  process  is  iterated  until  the  magnitudes  of  the  Zix^  are 
sufficiently  small  such  that  the  corrected  values  of  x^  are  said  to 
represent  the  navigation  solution. 


and 
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The  normal  equations  (6)  are  constrained  to  yield  a  solution  on 
the  surface  of  the  geoid  in  the  case  of  ship  navigation  or  on  the  surface 
of  a  co-geoid  in  the  case  of  aircraft  navigation.  The  constraint  is  applied 
by  means  of  the  theoretical  doppler  evaluation  on  the  geoid  or  co-geoid  while 
the  partial  derivatives  allow  only  movement  of  the  observer's  position  and 
velocity  on  the  geoid  or  co-geoid  for  ship  and  aircraft  navigation  respec¬ 
tively.  The  explicit  constraint  of  the  state  vector  r^j  is  given  by 
equation  (8). 

The  general  form  of  the  partial  derivatives  needed  for  the  normal 
equations  is  given  by  (7). 


i  =  1,  2,  3,  4 


The  partial  derivatives  may  be  obtained  in  a  straightforward  manner  beginning 
with  equations  (2),  (8)  and  (U). 

The  navigator's  position  on  the  geoid  in  earth  fixed  cylindrical 
coordinates,  is 


(8) 


p  P  ^ 

=  R^Cl  +  (1  -  f)  -tan  ®(t)]  +  (h  +  H)  cos  cpCx) 

tan  t(t)  +  (h  +  H)  sin  cp(t) 


where  and  f  are  the  equatorial  radius  and  flattening  of  the  ellipsoid 
respectively,  h  is  the  altitude  of  the  antenna  with  respect  to  the  geoid, 

H  is  the  altitude  of  the  geoid  with  respect  to  the  reference  ellipsoid  and 
is  usually  defined  as  the  geoidal  height.  Geoidal  height  may  be  computed 
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to  a  sufficiently  accirrate  approximation  as  the  difference  between  the  geoid 
potential  function  U,  and  the  reference  ellipsoid  potential  function  V, 
divided  by  the  local  gravity,  here  taken  as  the  partial  derivative  of  V  with 


respect  to  r. 


(9)  H  = 


N* 

r  00  N  n  ^ 

‘’ots  *^3  %  *  L  \  K  =“  ■"  ^ 

n=5  _ n=2  m=l _ J 


1  +  3C2°P2°+  SO^P,;- 


-  .  (u)  R  cos  cp)”^ 

0^0  ^  e  o 


where  p"*  are  spherical  harmonics;  C™,  S  are  the  coefficients  of  the 
harmonics  and  p  is  the  product  of  the  Newtonian  gravitation  constant  g, 
and  the  mass  of  the  earth,  M.  We  use  the  geopotential  function  given  by 
Guier  and  Newton,  [.8]  .  The  usual  definitions  of  cp  and  \  used  in  evaluating 
H  are  geocentric  latitude  and  longitude,  .respectively.  Negligible  errors 
are  obtained  by  using  geodetic  latitude  and  longitude.  Geoidal  height  is 
a  well  behaved  function,  changing  slowly  with  respect  to  latitude  and 
longitude  increments,  and  for  this  reason  it  will  only  be  necessary  to 
compute  H  once  per  navigation  fix  provided  the  original  estimates  of  (p  and 
X  are  correct  to  within  2  degrees. 

The  altitude  of  the  antenna  with  respect  to  the  mean  sea  level, 
the  geoid  in  the  case  of  a  ship,  is  treated  as  a  measured  system  constant. 
This  completes  the  specification  of  the  constraints  r^j  and  r^  for  the 
case  of  ship  navigation.  Aircraft  navigation  requires  external  measurement 
of  h  to  give  the  estimate  of  the  state  vector  rj^.  Sufficiently  fast  sampling 
rates  of  h  will  allow  an  estimate  of  r^j  to  be  given. 

The  effect  of  altitude,  h  or  H,  errors  in  the  computation  of  the 
navigator's  constrained  radius  is  shown  on  Figure  6.  The  residuals  of  this 
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effect  were  given  on  Figure  1  and  shown  to  be  an  odd  function.  Again,  the 
effect  of  the  odd  function  residual  can  be  misinterpreted  in  the  least  squares 
adjustment  as  either  a  longitude  error  or  a  velocity  north  error. 

Figure  6  shows  the  longitude  error  due  to  height  errors  of  40 
meters  and  80  meters.  The  relative  shape  of  the  curves  shows  this  effect 
to  be  non-linear.  Generally,  this  effect  is  negligible  for  ship  navigation 
since  the  capabilities  exist  to  compute  the  navigator's  radius  to  an 
accuracy  considerably  better  than  ij-0  meters.  Aircraft  navigation  will  require 
careful  measurements  and  modeling  of  the  least  squares  procedure.  The  effect 
of  altitude  error  on  the  velocity  north  is  negligible,  amounting  to  a 
maximum  value  of  about  0.05  knots  for  an  altitude  error  of  60  meters. 

The  observer's  inertial  longitude  is  given  by  (lO) 

(10)  A(X,t)  =  X(t)  +  +  Wg(T  -  t^) 

(t  ) 

where  Aq  °  is  defined  as  the  right  ascension  of  Greenwich  meridian  at 
an  epoch,  here  t^,  zero  hours  UT,  and  oo^,  the  earth's  average  sidereal 
rate. 

In  order  to  obtain  the  partial  derivatives,  we  decompose  equation 
(8)  into  inertial  cartesian  coordinates  using  (lO)  and  differentiate  as 
necessary.  The  partial  derivatives  needed  will  become  clear  if  the  first 
order  expression  for  the  doppler  shif't  is  expanded  as 

[xx  +  yy  +  zz]p  ^ 


CONFIDENTIAL 

Tlii>  document  containt  infonuatioa  aflfctihg  the  n%donal  defend  of  the  United  States  within  the  meanini;  of  the  Espionage  Laws,  Title  18 
U.S.C..  Sections  793  and  i94.  The  transmission  or  (he  revelation  of  ita  contents  in  any  manner  to  an  unaulhorited  person  is  prohibited  by  law. 


THE  HOPKINS  UNIVKKSITY 

APPLIED  PHYSICS  LABORATORY 

SlLVCM  SMIH4.  MAMYCANO 


-16- 


CONFIDENTIAL 


and  the  general  partial  derivative  is  then 


h. 

1 


-1 


±  =  1,  2,  3,  ^ 


The  least  squares  adjustment  of  the  data  is  normally  referenced 
to  an  epoch  since  in  the  general  case  the  latitude,  longitude,  and  velocity 
components  are  functions  of  bime,  t.  The  epoch  chosen  is  t^,  the  observer's 
time  of  satellite  closest  approach. 


(12) 


T 

<p(t)  =  cp(t^)  + 

T 

C 


T 

X(t)  =  \(t^)  +  I" 


q)(t)dt 


X(t  )dt 


Using  numerical  integration  of  oblate  earth  velocities  in  (32) 
permits  any  reasonable  maneuver  of  a  ship  or  aircraft,  depending  upon  the 
sample  rate  of  the  data.  The  maneuvering  capability  does  require  some  form 
of  an  external  dead  reckoning  device  such  as  an  inertial  system.  In  the 
event  no  maneuver  capability  is  required  during  the  time  of  satellite  pass, 
then  dead  reckoning  devices  are  not  necessary  for  furnishing  initial 
parameter  estimates. 

Equation  (6)  represents  a  positive  definite  quadratic  surface, 
as  is  evident  from  the  results  given  later.  The  quadratic  surface  is  highly 
elliptical  with  the  semi-major  axes  in  the  direction  of  longitude  and 
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velocity  east  with  the  minimum  of  the  surface  at  the  true  zero,  that  is 
where  errors  of  latitude,  longitude,  velocity  north  and  velocity  east  are 
zero.  Equation  (6)  is  positive  definite  when  the  doppler  data  used  in  the 
least  squares  adjustment  covers  a  sufficient  time,  or  geometrical  span. 
Figure  1  shows  the  latitude  and  velocity  east  residuals  to  be  virtually 
the  same  except  near  the  horizon.  Use  of  doppler  data  near  the  horizon 
introduces  another  systematic  error  due  to  tropospheric  refraction. 

Tropospheric  refraction  effects  can  contribute  an  error  which, 
in  general,  is  an  odd  function  with  r  espect  to  time  of  closest  approach, 

[9],  [lO].  Errors  behaving  as  odd  functions  will  be  misinterpreted  in  the 
least  squares  adjustment  as  false  longitude  or  velocity  north  residuals, 
resulting  in  errors  to  these  determinations.  Ordinarily,  tropospheric 
refraction  effects  need  not  be  accounted  for  if  the  portions  of  the 
doppler  curve  between  the  horizon  and  approximately  ten  degrees  elevation 
are  deleted  from  the  least  squares  computation.  However,  in  determination 
of  position  and  velocity,  it  has  been  found  necessary  to  use  the  doppler 
curve  at  elevation  angles  of  about  five  degrees.  Figure  1  shows  the 
doppler  residuals,  that  is  the  difference  between  the  theoretical,  or 
modeled,  doppler  data  and  the  measured  doppler  data  for  latitude,  longitude, 
velocity  north,  and  velocity  east  errors.  It  is  clear  from  Figure  1  that 
the  longitude  and  velocity  north  residuals  can  be  identified  as  odd  functions 
with  respect  to  satellite  closest  approach  while  latitude  and  velocity  east 
residuals  are  even  functions. 

The  velocity  error  effects  have  furnished  the  motivation  in 
seeking  to  determine  observer's  velocity  and  position  from  satellite  doppler 
data.  Until  recently,  only  partial  success  in  velocity  determination  had 
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been  achieved.  In  1^66  doppler  data  were  obtained  on  a  ship  maneuvering 
within  range  of  a  shore  based  theodolite  triangulation  net.  The  theodolite 
data  were  reduced  to  give  the  ship's  true  position  as  a  function  of  time. 

The  same  doppler  data  were  navigated  in  two  different  fitting  routines  as 

a)  Determination  of  frequency,  latitude,  and  longitude. 

b)  Determination  of  frequency,  latitude,  longitude,  and  velocity 
north. 

The  results  of  the  two  methods  of  navigating  are  shown  on  Figure 
7.  The  f^,  cp  and  X  determination  used  velocity  components  from  an  inertial 
navigation  system.  The  f^,  cp,  X  and  Vj^  determination  used  the  velocity 
north  component  from  the  inertial  system  only  as  an  initial  estimate  in  the 
least  squares  adjustment.  The  resultant  improvements  in  the  f^,  cp,  \, 
navigation  fixes  are  due  to  the  determination  of  a  more  accurate  velocity 
north  component  in  the  fitting  process  using  satellite  doppler  data. 

Recent  experimentation  by  G.  C.  Weiffenbach  and  the  author  indicates 
that  velocity  east  may  also  be  determined  provided  that  the  satellite 
oscillator  and  receiver  reference  oscillator  are  sufficiently  stable  such 
that  the  frequency  term  does  not  have  to  be  included  in  the  normal  equations 
as  indicated  by  equations  (3)  thru  (6).  Weiffenbach  and  the  author  have 
used  real  world  data  in  arriving  at  the  required  oscillator  stability.  The 
results  from  real  data  and  the  simulation  are  in  substantial  agreement. 

The  experiment  using  the  simulation  program  consisted  of  implementing 
equation  (6)  and  the  constraints  noted  along  with  controlled  errors  repre¬ 
senting  oscillator  instabilities.  Figure  8  shows  the  results  of  the 
simulation  experiment,  that  is,  the  effect  of  velocity  east  errors  due  to 
oscillator  instabilities.  The  results  shown  on  Figure  6  can  be  expressed  as 

|5Vj,|  =  |5f^l  8/Ea(*c)) 
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Klevatlon  Angle  at  Closest  Approach,  E,^(tQ),  degrees 

Fig.  8  Velocity  East  Error  per  Frequency  Bias  Error 
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A  value  for  the  relative  oscillator  stability  can  be  determined 
by  assigning  a  prescribed  value  to  the  allowable  velocity  east  error.  It 
is  more  meaningful  to  include  the  effect  of  &V„  on  the  latitude  error  using 
the  curve  of  Figure  4  as 

|8<p|  . 

Substituting, 

Istpl  =  |5f^|g^(Ejt^))fj,(Ejt^)) 

a  value  of  |8f^(/:^,  selected  to  allow  0.01  nautical  mile  latitude  error, 

12 

represents  a  total  oscillator  instability  of  5  parts  in  10  .  This  total 

oscillator  instability  includes  the  combined  effects  of  satellite  oscil¬ 
lator  and  receiver  reference  oscillator  instabilities. 

Similar  results  have  been  obtained  with  real  data  in  the 
experiment  by  Weiffenbach  and  the  author  using  navigation  results  from  two 
separate  fixed  sites.  One  fixed  site  is  equipped  with  a  Cesium  frequency 
standard  while  the  other  uses  a  Rubidium  reference.  At  the  time  of  this 
writing,  the  Rubidium  standard  appears  superior  in  short  terra  stability 
over  a  20  minute  span. 

Assuming  that  it  is  possible  to  acquire  atomic  frequency  standards 

12 

with  a  stability  of  about  2  parts  in  10  for  both  the  satellite  and 
receiver  reference  oscillators  then  it  is  probable  that  velocity  east  can 
be  determined  in  the  real  world  with  good  accuracy.  Determination  of 
velocity  north  has  already  been  demonstrated  on  an  experimental  basis, 
however  it  must  be  noted  that  acquisition  of  a  value  for  velocity  north 
places  rather  stringent  requirements  on  satellite  velocity.  Equation  (2) 
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clearly  states  that  the  ephemeris  is  the  navigation  reference.  Satellite 
velocity  must  he  accurate  to  less  than  0.05  meters  second  if  a  0.1  knot 
determination  of  velocity  north  is  desired  since  one  knot  is  0.51*<^  meters 
per  second. 

III.  CONCLUSION 

Position  and  velocity  have  been  determined  in  a  computer  simulation. 
The  simulation  program  accounts  for  all  modeling  except  atmospheric,  solar, 
and  magnetic  activity.  Excluded  are  ionospheric  refraction  [l],  [3], 
tropospheric  refraction  C9]>Cl0]>  and  the  effects  of  drag  forces  which  must 
be  accounted  for  in  the  real  world.  Assuming  sufficient  oscillator 
stability  and  a  sufficiently  accurate  ephemeris,  then  the  accuracy  of  the 
doppler  satellite  navigation  system  may  be  theoretically  indicated  by  the 
computer  simulation  results  given  on  Figure  9* 
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Orbit  Velocity  Error  <  0.05  meters/sec 
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Oscillator  Stability  2(l0) 
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Fig.  9  Theoretical  Navigation  Accuracy 
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